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ABSTRACT 

A generalized mean magnetic field induction equation for differential rotators is derived, 
including a compressibility, and the anisotropy induced on the turbulent quantities from the 
mean magnetic field itself and a mean velocity shear. Derivations of the mean field equations 
often do not emphasize that there must be anisotropy and inhomogeneity in the turbulence 
for mean field growth. The anisotropy from shear is the source of a term involving the 
product of the mean velocity gradient and the cross-helicity correlation of the isotropic parts 
of the fluctuating velocity and magnetic field, (v • b)(°). The full mean field equations are 
derived to linear order in mean fields, but it is also shown that the cross-helicity term survives 
to all orders in the velocity shear. This cross-helicity term can obviate the need for a pre- 
existing seed mean magnetic field for mean field growth: though a fluctuating seed field is 
necessary for a non-vanishing cross-helicity, the term can produce linear (in time) mean field 
growth of the toroidal field from zero mean field. After one vertical diffusion time, the cross- 
helicity term becomes sub-dominant and dynamo exponential amplification/sustenance of 
the mean field can subsequently ensue. The cross-helicity term should produce odd symmetry 
in the mean magnetic field, in contrast to the usually favored even modes of the dynamo 
amplification in sheared discs. This may be important for the observed mean field geometries 
of spiral galaxies. The strength of the mean seed field provided by the cross- helicity depends 
linearly on the magnitude of the cross-helicity. 

Subject Headings: magnetic fields; galaxies: magnetic fields; turbulence; accretion 
discs. 
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1. Introduction 

Mean magnetic field dynamo theory has been a leading formalism to characterize the 
evolution and origin of large-scale magnetic fields in stars and Galaxies (e.g. Moffatt 1978; 
Parker 1979; Krause & Radler 1980; Zeldovich et al. 1983). The mean field dynamo appeals 
to a combination of helical turbulence, differential rotation, and turbulent diffusion to ex- 
ponentiate an initial seed mean magnetic field. The standard textbook kinematic theory is 
widely known to be incomplete, because it ignores the back reaction of the growing magnetic 
field on the turbulence driving the field growth. 

The extent of dynamo growth quenching from the backreaction is not agreed upon 
(Cowling 1957; Piddington 1981; Zeldovich et al. 1983; Vainshtein & Cattaneo 1992; Catta- 
neo & Hughes 1996; Vainshtein 1998; Field et al. 1999). Field et al. (1999) and Blackman & 
Field (1999ab) have suggested that some previous criticisms of enhanced suppression might 
be challenged. Blackman & Field 1999b argue that simulations which employ periodic bound- 
ary conditions (e.g. Cattaneo & Hughes 1996) cannot test for suppression of dynamo theory 
because an upper limit on the required dynamo quantities can be shown to be strongly 
restricted when boundary terms are ignored. While some simulations done with diffusive 
boundary conditions do show evidence for a mean-field dynamo in fully non-linear MHD 
(Brandenburg & Donner 1997), there also exist non-linear analytic models (Vainshtein 1998) 
which show more extreme suppression. The restrictions with respect to actual astrophysical 
systems needs to be investigated and more analytic and numerical studies are pending. 

Since the backreaction consequences are not fully resolved, and because some results 
of backreaction studies do not support catastrophic quenching of the dynamo theory, there 
remains ample motivation to continue studying the solutions of the mean-field dynamo equa- 
tions in astrophysical systems. An industry continues to do so, motivated by some success 
in the solar (c.f. Parker 1979; Belvedere 1990), galactic (c.f. Ruzmaikin et al 1988; Beck et 
al 1996), and accretion disc (Brandenburg & Donner 1997) cases and the simplicity of the 
theory (c.f. Moffatt 1978; Parker 1979). The focus is often on the solutions, rather than the 
derivation of the equations being solved. 

In this paper, I derive a generalized dynamo equation, including a restricted compress- 
ibility and the anisotropy induced from the backreaction of the mean magnetic field and from 
a mean velocity shear. The inclusion of the shear induced anisotropy leads to an important 
term in the mean field evolution equations, involving the product of mean velocity gradients 
and the cross-helicity. This term should be included in dynamo models of sheared systems 
such as galaxies, stars, and accretion discs, as it turns out to play a fundamental role in 
the mean field growth as long as the cross-helicity is non-vanishing. In particular, it can 
generate a mean magnetic field when the initial mean field is zero as first pointed out by 
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Blackman & Chou (1997) and Brandenburg & Urpin (1998). 

Previous studies have recognized the potential importance of the cross-helicity term. 
Yoshizawa and Yokoi (1993) and Yokoi (1996) first discussed the term in the astrophysical 
context. Brandenburg & Urpin (1998) have produced a nice analysis of the dynamo equation 
with shear. The results herein are closer to those of Brandenburg & Urpin (1998). However, 
there is an ambiguity in their derivation of the employed dynamo equations of the type 
emphasized by Blackman & Field (1999a). It is important to distinguish the isotropic com- 
ponent of the turbulence from the anisotropic component, both of which must be present. 
The formalism herein explicitly avoids this ambiguity and provides a more complete deriva- 
tion of the generalized dynamo equations with the important shear term. In addition, it will 
be shown that the cross-helicity term survives to all orders in the shear not just to linear 
order. 

In section 2 the basic mean field equations are derived including anisotropics induced 
both from the mean magnetic field and the velocity shear and compressibility. In section 3, I 
show that the term which will lead to linear mean field growth from zero mean field survives 
all orders in the shear anisotropy. In section 4 the mean magnetic field equation is solved in 
the early time limit when the shear term dominates. After one vertical diffusion time from 
the chosen t — 0, this term becomes sub-dominant. The resulting solutions are then given 
for both a — Q and a 2 dynamos. Section 5 discusses the implications of these solutions, with 
emphasis on galaxies, and section 6 is the conclusion. 

2. Mean field evolution equation 

Mean field dynamo theory characterizes the growth of magnetic fields on scales larger 
than the input turbulence. The induction equation describing the evolution of the magnetic 
field is given by 

d t B = V x (V x B) + AV 2 B, (1) 

where A is the micro-physical diffusivity and the pressure gradient and is taken parallel 
to the density gradient to ignore the Biermann-Battery type term (Biermann 1950). The 
Navier-Stokes equation describing the velocity evolution is given by 

pd t V = -V • VV - VP - V(5 2 /8tt) + B • VB/(4tt) + F(x, t) + V0, (2) 

where V is the velocity, B is the magnetic field, P is the pressure, and F is a forcing 
function, and V0 is a gravitation term. Splitting the magnetic field into mean and fluctuating 
components we write B = b+B where (b) = 0, and the velocity V is similarly defined. (Note 
that since astrophysical rotators usually have turbulent input scales that are a non-trivial 
fraction of the mean scales, it is important to distinguish between fluctuating quantities and 
small scale quantities. Their magnitudes differ by a quantity that varies as the 3/2 power of 
the ratio of the turbulent input scale to the mean scale.) 
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We assume that the mean velocity is time independent, so we do not consider the mean 
velocity equation any further. For the mean induction equation we have (e.g. Moffatt 1978) 

d t B = V x (v x b) + V x (V x B), (3) 

where the brackets and the over-bar indicate mean values and the lower case indicates fluc- 
tuating values. 

To write the fluctuating components, I follow the approach of Field et al. (1999) and 
Blackman & Chou (1997). The fluctuating fields are written in terms of the zeroth order 
isotropic base state and higher order anisotropic components which are due to the presence 
of a mean field. Here we consider the influence of both the mean magnetic field as well as the 
mean velocity field. This approach can be motivated when the magnitude of the fluctuating 
anisotropic component of the velocity and magnetic field are < the mean velocity and mean 
magnetic field respectively. 

The zeroth order quantities are not solved for but have the properties of isotropic ho- 
mogeneous turbulence. The base state is also taken to have a reflection asymmetry. In 
astrophysical systems this property can result from the presence of an underlying rotation 
and a density gradient (Moffatt 1978; Krause & Radler 1980). One can take the liberty of 
presuming that the forcing function acts only on the zeroth order state and feeds it with 
isotropic but reflection asymmetric turbulence. Since we do not solve the zeroth order equa- 
tions, we do not need to specify the explicit form of the forcing function. The anisotropy 
induced from the mean velocity and mean magnetic field will be considered in the higher 
order components. Compressibility is allowed for in a restricted sense: the density will be 
taken to have a spatially dependent mean component p, and a fluctuating component which 
is free of influence of the mean fields. 

Subtracting @ from (|l|) we have for the fluctuating field 

d t b = V x (v x b) - V x (v x b) + V x (v x B) + V x (V x b) (4) 
For the fluctuating Navier-Stokes equation we have 

pdtv = -v ■ VV + (v • Vv) - v • Vv - Vp - V(B • b)/47r - W/8tt + V(6 2 )/8tt 

+F(x, t) + B • Vb/4vr + b • VB/4tt + b • Vb/47r - (b ■ Vb) /4vr. ' ° ' 

We have ignored the microphsical viscosities for present purposes, as they are taken to be 
small. Write b = b^ ^ + and similarly v = v^ ) + v^ A ' where and are the 
anisotropic parts of the fluctuating components. As long as |b( A )|(~ 47rp|v( A )|) < |B| and 
|V|, then v^ A ) and b*^ can in principle, be solved for explicitly in terms of all orders in V 
and B (see Field et al. 1999 for the solution to all orders in B ignoring mean field gradients). 
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In this section I restrict the calculation to linear order in both |B|/|b| and (|V|/i?)/(|b|/L), 
which is sufficient for the main points, and then in the next section I expand to all order in 
(|V|/i?)/(|b|/L), since this quantity is of order 1 for many applications. 

Working in the local standard of rest (LSR) frame (Field et al. 1999) in which V = 0, 
but not VV to first order in the mean fields we have 

O t bW = b(°) ■ VvW + b« ■ Vv<°) - b(°)V ■ v« - b«V • v(°) - V x (v x b)W 

_ v (0) . Vb (l) _ v (l) . Vb (0) + v (0)y . b (l) + v (l)y . b (0) (g) 

+ B ■ Vv< ' - BV • v(°) - v(°) • VB + b< > • VV - b^V • V, 



and 



pd t v« = -v« • VV + (v ■ Vv)« - yM • Vv« - v« ■ VvW 
_VpW - v (°) ■ VvW - J B i VbS 0) /47r - 6$V&«/47T - 6«V6^/4tt 
-b$VB m /4ir + V(6 2 )( 1 )/4tt + B • Vb(°)/47r + b<°> • VB/4tt+ 

(b ■ Vb)W/47T - (b ■ Vb)«/47r. 



(7) 



Now 



(v x b)W = ( v (°) x J d t h^dt) - (bW x J d t w^dt). (8) 
Using (|^) and (|7|) in (§) we have for the first term on the right of ([8]) 

(v<°> x Jdtb^dt), = r c {e ijk vfd m v^)B m - T c {e ijk vfd m v$)B k 

-T c (e i]k vfv^)d m B k - r c (e i]k vfd^b k 0) )d m V m (9) 

+T c (e ijk vf ) b$}d m V k - Tc^ijkvfh^dmVm- 



After isotropizing (isotropizing means that 2 and 3-index tensor average correlations of ze- 
roth order quantities are proportional to 5y and tensors respectively) and making the 
controversial (e.g. Cattaneo 1994) but commonly employed replacement the time integrals 
with multiples of the correlation time, we have 



,(0) 



/ d t h^dt) = -((v(°' • b(°>)V x V - (v(°> • v<°))V x B - (v<® ■ V x v<°>)B) (10) 
J 3 



3 

Following the same procedure as above, I obtain for the second term on the right of 
-(b(°> x fdtv^dt) = f ((bW • v(°))V x V + 2(b^ • V x b(°))B/(47rp) 



(11) 



_ 2 (b(°) ■ b(°))V x B/(47rp) + p-^b x Vp (1) )). 
Expand the last term in fllTD to first order: the energy equation can be written 

d t P + V • VP = 7PV -v + C, (12) 
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where 7 is the adiabatic index and C is a cooling or heating term which we assume to have 
only a mean contribution. Writing the fluctuating part of this equation to first order in the 
mean fields then gives 

%,« = - v (0) • VP+ (v ■ Vp) (1) - ( v (0) • Vp) (1) +7PV ■ v (0) + 7 pV • V +7(pV • v) (1) -7 (pV ■ v) (1) . 

(13) 

Plugging this back into the last term of (|Tl~D and using V ■ V = gives 



p-^b x Vp^), = {(e ijk bfd k Jd tP Vdt>) ~ r c ((-ey fc (6f ^9^°))^, - e ijk (bf d k v$)d m P 

-e ijk (bfv$)d k d m P 
+je ijk (bf ) d k d m v$)P + je^ibfd^dkP). 

(14) 

The 1st 3rd and 4th terms on the right of ( P) ) vanish from anti-symmetrization after 
isotropization. The 5th term vanishes directly from isotropy of zeroth order correlations. 
Thus the only remaining term is the 2nd term. Combining this term with (ID) and ([TTJ) 
gives 

(v x b)W = g(b • V x v)(°»VP + (r c B/3)(2(b • V x b)(°V(4vrp) 
-(v V x v)W) - r c V x B((b • b)(°7(67rp) + f(v ■ v)^) + x V((v • b)(°>) (15) 

ee «S 0) VP + c4 0) B - {3{ 0) V x B + V x V 

Taking the curl of (|15"D, plugging into (0), and again assuming Vp x VP = the equation 
for the evolution of the mean magnetic field becomes 



d t B = T c ((b ■ V x b)(°V(67rp) - (v ■ V x v) (0) /3)V x B 
+r c ((b ■ b)(°)/(67rp) - (v ■ V x v)(°)/3)V 2 B - T c /(Qnp 2 )(b ■ V x b)^Vp x B 
+r c (b ■ b)(°>Vp x(VxB) + |(b. v)<°'V xVxV + Vx(VxB) 
= 4 0) V x B + Vaf x B + /3{ 0) V 2 B - V^ (0) x (V x B) - ^ 0) V 2 V + B-VV 



(16) 



The last term of ( |T6"D is the standard f2 shearing term and the c4 and /S^ '' terms are the 
standard a and /3 dynamo coefficients (Moffatt 1978) except for the inclusion of a spatially 
dependent density. We will see that the penultimate term on the right can generate a seed 
field. The coupling coefficient (3^ is the cross-helicity. The importance of a cross-helicity 
term was first addressed by Yoshizawa and Yokoi (1993), Yokoi (1996) but not as a source 
of seed field for the standard dynamo. Blackman & Chou (1997) and Brandenburg & Urpin 
(1998) first pointed out the seed field role. Before solving the generalized mean- field equation 
which includes the role of the cross-helicity, I will show that the cross-helicity term survives 
to all orders in the shear. 

3. Survival of the linear growth term to all orders in mean shear 
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The expansion in section 2, was taken to linear order in |B|/|b (0) | and (|V|/i?)/(|b (0) |/L), 
where R is the scale of variation of V and L is the turbulent outer scale. While | B | / 1 b^ ^ | < 1, 
<< 1 for early dynamo evolution, and < 1 in the Galaxy at present, and arguably < 1 in 
accretion discs, the quantity (| V|/i?)/(|b^|/L) ~ 1 in accretion discs for when a shearing 
instability drives the turbulence (Balbus & Hawley 1998). This is simply the statement that 
the instability growth time is of order the rotation time. For this purpose it is necessary to 
expand to all orders in at least (\\\/R)/(\b^\/L). 

In this section I show that the second last term found in ( ]T6|) survives to all orders in the 
mean shear field when the V is dominated by 1 component (V^ in the present case). Splitting 
up the small scale fields into isotropic and anisotropic components we recall b = b^ ) + 
and v = v 1 - - 1 + v^. We then have for the turbulent EMF 

(v x b) = (v<°> x b^> + (v^ x b< >> + (v< A > x b^>. (17) 

To simply show that the linear growth term survives to non-linear order in the shear, let us 
take the mean magnetic field to be identically zero and ignore the pressure term. (This will 
have to be generalized in the future, but the issue of survival of the term in question is not 
affected.) Then, the evolution equations for the nth order anisotropic components become 
(ignoring second order correlations in r c ) simply 

cW n) = -v^djVi, (18) 

and 

= ///' '^r, (19) 

Summing over n we have 

d t v\ A) + vfdjVi = -vfdfr, (20) 

and 

dtb^ -bf ) d j V i = bfd j V i . (21) 

Multiplying by an integrating factor, and defining M iq {t) = [e vv *]j 9 (where VV is indepen- 
dent of time) these equations can be written 

dtiv^M^t)) = -vfdjViMit)^, (22) 

and 

dtib^ M iq {-t)) = bfd^M^-t). (23) 
The solutions to these equations are given by 

= v^M^M-^t) = - J* vfd^M^dt'M-^t) + v s (0), (24) 
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and t 

b[ A) = M iq {-t)M- s \-t) = j o bfd&Mni-t'WM-^-t) + 6.(0) (25) 

The 3 terms in the turbulent EMF ([T7|) then become 

(v (0) xb^} p = (e pms v^ J &f M*(-OdO0;WCH) = MS^f)^ = |(v (0) -b^>(Vx V) p 

(26) 

and 

-(b(°)xv^) p = (e pms bM J vfM^tyt'^M-^t) = (e^v^bf^V, = ^(^.b^)(Vx V) p 

(27) 

and finally 
(v (A) x b^> p = (<^(- jjvfajt^ 

(e pms ^(0)M im (t)6; A) (0)M JS (t)). 

(28) 

Notice that the sum of (p6|) and ( p7|) is simply /Jj which enters in (|T5l) and thus (|T6|) . 
Thus ( p8|) appears to produce an extra term in the EMF. However, it vanishes when V is 
dominated by 1 component (in our case Vj ~ Vr, the toroidal mean velocity): the index j 
in Mij(t) always corresponds to the index of Vj when the exponential in My is expanded. 



This means that m = s in (p8[) and due to the e mps , both terms in (28) vanish. Thus, we see 
that in the case for which the mean velocity is dominated by one component, the shearing 
term contributes a term which will lead to linear growth from zero mean field, as per the 
next section, even when the anisotropy due to this shear is included to all orders. 

3. Growth of the Mean Field 

For simplicity, I will assume that the mean field is axisymmetric. All of the components 
of the magnetic field can then be expressed in terms of the toroidal field and the poloidal 
field V x A?, where A is the vector potential. The equation for the poloidal mean field can 
then be "uncurled" (since an arbitrary gradient term vanishes for axi-symmetry) . From ([16]) 
we then have 

d t A T = c4 0) B T - /?i 0) V x V x A T + ^ 0) (V x V) T + (V x B) T . (29) 
For the toroidal component of B we have 

d t B T = -c4 0) V 2 A T - (Vc4 0) -) A T + (3[ 0) V 2 B t + (V x A t ) • VV + (V(3[ 0) ■ V)B T - ^ 0) V 2 V T . 

(30) 

Let us now proceed to solve equations (EH) and (|3"0"D. 



3a. Linear Growth Regime 
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At early times, the dominant term on the right of (^) is the term Q(r, z) = (3^ V 2 V, 
taken to be independent of time. Solving (j30D with only this term gives 



B T = -Q(r,*)t, (31) 

if we assume the initial B is negligible. Plugging this solution into (|29|) , and solving we 
obtain, to lowest order in t 

A T ~ 2H 2 a { 2 ] Q(r } z)t/(3{ 0) . (32) 

where H is the characteristic gradient length of the mean magnetic field. How long does this 
solution hold approximately? A typical value of the other terms in (p9[ ) and (|3~0D is given 
by the diffusion terms. Thus plugging ([H]) or (|32"D into the respective diffusion terms, with 
the strongest variation of the magnetic field taken to be the vertical direction, one finds that 
the linear growth term becomes subdominant at t > H 2 /ft[°\ At this time, the fields have 
grown to the values 

»•/,,-./ ^ -(#7/3i 0) )4 0) V 2 V T (33) 

and 

A T , mt ~ H*/(pF>)*a£> /3f V 2 V T . (34) 

Because the fluctuating magnetic field builds up to near equipartition with the turbulent 
kinetic motions with a growth rate ~ eddy turnover time, |b^|/(47rp) 1 / 2 ~ |v(°)| at times 
of interest for long term mean field evolution. We then have 

V A ,cru~ \V T \(H 2 /R 2 )x i0 \ (35) 

where VA,crit is the Alfven speed associated with the mean magnetic field, — [|/^2°' ) |(4vrp) _1 / 2 ]/|/3 1 ' ) | 
and VA,crit is the Alfven speed associated with the mean magnetic field at the time when 
the linear term becomes subdominant. Note the importance of the ratio of cross-helicity 
magnitude to turbulent diffusivity magnitude. 

3b. Dynamo Growth Regime 

Now I proceed to solve fl29| ) and ( p0|) in the in the t > H 2 / limit, when the when the 
third term on the right of (|2"P|) and the last term on the right of (BUI) can be ignored. In this 



regime, writing (|29D and (|30| ) in terms of their scalar components in cylindrical coordinates 
gives 

dtA* = o4 o) S + P{ 0) V 2 A^ - (3[ 0) A^/r 2 , (36) 
where indicates the toroidal component, and 

d t B+ = -o$ ) ^A++a$ ) A+/i*-Vo§ ) -VAi+^^^^ 

(37) 
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I assume the case where the vertical gradients of B dominate and the radial gradients of 
Vfj, dominate and look for solutions of the form cx exp[nt+ik-x] and, =oc exp[rji+ik-x]. 
(Given the spatial density gradients in a system, one can be more accurate, since the A^^ rit 
and B^u have spatial dependences which depend only on the density gradient. For present 
purposes I take the simpler approach.) From (|36D and fl37j) , we obtain 

(n + /3{% 2 )^-4% = 0, (38) 

and 

(a^k 2 - ik z d z af - ik z d r %)A^ - (n + p{ 0) k 2 - ik z d z p[ 0) )B^ = 0. (39) 
The dispersion relation resulting from (|38|) and ( |3~9"D is 

n 2 + n(2(3[ 0) k 2 z - ik z d^) + {K^f ~ ~ (*,a<°>) a + iafk z d z af + % afk z d T %. 

(40) 

The solution to this quadratic is 

2n = -k z (3[ 0) + ihjdtffl ± [Ak 2 z af 2 - k 2 z {d z ^f - iAk z o$\d x af> + d^)} 1 ' 2 . (41) 

Scaling the terms inside the square brackets on the right of fl4"T|) shows that the second term is 
down from the first by some constant < 1 times (L/H) 2 , where L is the dominant turbulent 
scale. We assume that this second term can be ignored for simplicity. The remaining real 
part of fl4l"l) is of interest. This can be found from the following: write the right side as 
(c + di) 1 / 2 , where c and d are real and seek the real quantity a such that a + bi = (c + d%) x l 2 . 
In general, this implies that 

2a 2 ={c±{c 2 + d 2 ) 1/2 }. (42) 

Applying (g^) to (gT]) gives 

Re[n] = -k 2 p{ 0) ± ^(4 0) /2^)[i ± [i + + d r %f / [Ka^f] 1 ' 2 ] 1 ' 2 . (43) 

There are two interesting limits of fl43|). In the limit that the differential rotation 
dominates the inner sum, we have 

Re[n] ~ -k 2 j{ 0) ± (fc,4 0) ) 1 / 2 (9 r ^) 1 / 2 /2 1 / 2 . (44) 

This has growing modes for k z < (\a^ d r V<p\/2l3^ ) 1 / 3 . The maximum growth wavenumber 

is k z , max = (|4 0) d r ^2 5 /^ 0)2 ) 1/3 so that n max ~ 0.3(|a^ 0) 9 r ^| 2 M 0) ) 1/3 . This is an "a-fi" 
type dynamo (c.f. Moffatt 1978; Ruzmaikin et al. 1988). 

In the limit that d z a2^ dominates or is comparable to <9 r V^, we have 
Re[n] * -A; 2 M 0) ±^4 0) [l±[l+(^4 0) /^4 0) ) 2 ] 1/2 ] 1/2 -k 2 z (3 { ?±k z af[l±2-" 2 f 2 , (45) 
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where the latter similarity follows when \(d z a 2 )/(k z a 2 )\ ~ 1. This has growing modes 
for k z < l.3a ( 2 0) /p{ 0) . The maximum growth wave number is at ~ 0.65a^/f3[°^ so that 
n max ~ 0.4(4 0) ) 2 //?f } . This is an "a 2 " type dynamo (c.f. Moffatt 1978). 

4. Discussion 

The results of sections 3a and 3b show the two phases of mean field growth for a shearing 
(disc) rotator. The first regime does not require a seed mean field, only a fluctuating field 
whose mean can vanish. The toroidal field produced in the first phase (the linear growth 
regime) should be asymmetric with respect to the disc plane (like an AO mode, e.g. Beck 
et al. 1996) since the cross- helicity (3 { 2 ] = (b(°) • v(°)) is a pseudo-scalar. The poloidal field 
growth feeds on this linear toroidal field. If the cross-helicity has the same sign throughout 
a hemisphere, the linear growth phase should not produce radial field reversals in that 
hemisphere since the mean field vector is determined only by the sign of the cross-helicity 
and the direction of the mean velocity and its gradient. 

The initial growth phase is important since standard dynamo models generally consider 
the equations only in the second phase. In the second phase it is generally easier to amplify SO 
type modes, whose toroidal field is symmetric with respect to the galactic plane (Ruzmaikin 
et al 1988; Beck et al 1996), though more complicated models are able to excite mixed modes 
e.g. AO-SO (Moss & Tuominen 1990; Moss et al. 1993). The dominant observed field after 
amplification saturates also depends on which seed field modes were present to begin with 
(Ruzmaikin et al. 1988; Poezd et al. 1993). Suppose an asymmetric seed field is present 
with magnitude well in excess of the seed symmetric mode. Then even if the asymmetric 
mode grows more slowly than the symmetric mode, the difference of initial magnitudes can 
more than compensate and the asymmetric mode could dominate. Also, once equipartition 
is reached, different modes may oscillate with different frequencies. 

In the Galaxy for example, the magnitude of the asymmetric seed resulting from the 
linear growth phase is, assuming no vertical shear and a radius of 8Kpc, from (^) B ~ 
2.7 x 10~ 9 (x (0) /0-01)( J ff/0.5Kpc) 2 ( J R/8Kpc)- 2 (^/2 x 10 7 cm/s)G. This is a substantial seed 
field, which would arise in one vertical diffusion time. I have scaled with x (0) ~ 0.01. 
However, note in the Galaxy for example, the dynamo can exponentiate > 20 or so times 
in the Galactic lifetime. Thus even if x ~ 10 -5 the seed from the linear growth phase 
could be important for seeding the subsequent dynamo growth up to its present value of few 
xl0~ 6 G. If — 0.01, or if there were vertical shear, then fl3~5|) would imply an even larger 
seed, and a much shorter time for the cross-helicity + dynamo to produce an equipartition 
mean field. The exponential amplification rate of a symmetric mode from a seed field of 
10~ 9 G would have to be larger than that of the asymmetric mode to be competitive with 
that seeded from the asymmetric linear growth phase. 
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Han et al. (1998) suggest that the Faraday rotation data for M31 supports the presence 
of an even mode (SO). For our Milky way however, the data favor a dominant AO mode (Han 
et al. 1997). While more galaxy data are needed, there appears to be the possibility that 
different modes dominate in different galaxies. The results herein suggest that the cross- 
helicity term can be important in determining the field geometry in the framework of in situ 
field generation models. 

A number of simple but important complications affect the observational implications of 
the above results. First, inhomogeneity and large scale local structures can lead to coherent 
local structures. There may also be seed fields present from protogalactic, cosmological, or 
supernovae origins which are amplified concurrently. In addition, the mean field in a galaxy 
is likely never to be zero since the dominant input scale of the turbulence in galaxies is always 
a non-trivial fraction of the galactic radius. For example, at the solar location of our Galaxy, 
the scale ratio is at most (lOOpc/lOkpc). Thus, just the random RMS field (Blackman 1998), 
or RMS contributions to dynamo coefficients that result (Vishniac & Brandenburg 1997) may 
be important in influencing the observation of a mean field. (It is important to note that 
mean field theory and Faraday rotation measurements are degenerate with respect to mean 
field topology: they cannot distinguish between a large scale field formed from averaging 
over small scale loops, or from that formed by averaging of a topologically connected field 
line with fluctuations. Spectral energy distribution approaches also share this degeneracy.) 

The two phase growth described herein, as that of Brandenburg & Urpin (1998), could 
be testable in some non-linear MHD disc simulations. Note that periodic boundary condi- 
tion simulations cannot test for mean field growth because the mean field is conserved by 
construction (c.f. Balbus & Hawley 1998). It is known that diffusion through the boundary 
is required for a working dynamo (Parker 1979; Ruzmaikin et al 1988), and it has also been 
realized that diffusive boundary conditions are required for a non-vanishing effect even 
when dynamo action is not present (Blackman & Field 1999b). There is some evidence that 
a dynamo is operating in steady accretion disc simulations with the appropriate diffusive 
boundary conditions (Brandenburg & Donner 1997). 

Brandenburg & Urpin (1998) point out that some previous simulations have estimated 
X- (Note that this is not x^ which comes in directly into the present formalism, but the 
full X-, to all orders in the mean fields. The two should be distinguished.) It is found that 
X ~ 0.03 in stratified convection simulations (Brandenburg et. al. 1996), x ~ 3 x 10~ 4 
in magneto-shearing driven instability simulations (Brandenburg et. al. 1995), and x ~ 
5 x 10~ 3 in supernova induced turbulent flows (Korpi et al. 1998; though these have not 
reached the saturated steady-state) potentially relevant for galaxies. Chandran & Rodriguez 
(1997) studied the evolution of the cross-helicity in the framework of the direct interaction 
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approximation, and do not rule out the possibility of significant cross-helicity on the turbulent 
input scale. Future simulations should give a better handle on \ an d X^°\ an d a comparison 
of the two in the various settings. As pointed out earlier, even values of 10~ 5 might be 
significant. Note that it is also important to distinguish between an RMS value that may 
fluctuate in time and a systematic value that remains steady in time. 

The use of a linear backreaction model is certainly incomplete (as are most studies of 
dynamo equation solutions). Both the shear and B should be included to all orders. Field 
et al. (1999) solved for all orders in B without field gradients. Ultimately, more analytic 
and numerical work are needed to understand how well the dynamo survives in the non- 
linear backreaction. Interestingly, dynamo theory in its standard form is the "complement" 
of analytic accretion on disc theory: the latter ignores the dynamics of the magnetic field, 
while the former does not fully include the magnetic backreaction on the velocity dynamics. 
Both require some closure approximation to turbulence to make analytic progress. Studies 
which focus on the solutions of the accretion disc or mean field dynamo equations often gloss 
over the approximations which led to the form of the equations being solved. 

5. Conclusions 

The equations for the evolution of the mean magnetic field in a sheared rotator have 
been derived, taking into account the anisotropy induced in the turbulence from the large 
scale magnetic field and the differential shear. A restricted compressibility has also been 
included. A shear term violating the homogeneity of the dynamo equations was shown to 
survive to all orders in the anisotropy. For finite cross-helicity, the mean magnetic field 
can grow in two phases. Linear growth ensues for one vertical diffusion time. The mean 
toroidal and poloidal fields grow without a seed mean field to a value whose Alfven speed 
is of order > (H/R) 2 x^V^ in an AO mode (toroidal field anti-symmetric with respect to 
the disc plane). After the vertical diffusion time, the second phase incurs, and the more 
"standard" dynamo evolution takes over. If the linear growth phase seed field dominated 
the field geometry, there would be a reversal in the toroidal field across the disc plane and no 
reversals in radius if the cross-helicity has the same sign in a given hemisphere. Determining 
the magnitude of x^ an d the full x to all orders for various applications will continue to 
be aided by numerical simulations. For non-vanishing cross-helicity, equation (16) should 
replace the standard textbook dynamo equations for sheared rotators. 
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